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The Poynting vector is an invaluable tool for analysing electromagnetic problems. However, even 
a rigorous stress-energy tensor approach can still leave us with the question: is it best defined as 
E X H or sls D X Bl Typical electromagnetic treatments provide yet another perspective: they 
regard E x B as the appropriate definition, because E and B are taken to be the fundamental 
electromagnetic fields. The astute reader will even notice the fourth possible combination of fields: 
i.e. D X H . Faced with this diverse selection, we have decided to treat each possible flux vector 
on its merits, deriving its associated energy continuity equation but applying minimal restrictions 
to the allowed host media. We then discuss each form, and how it represents the response of the 
medium. Finally, we derive a propagation equation for each flux vector using a directional flelds 
approach; a useful result which enables further interpretation of each flux and its interaction with 
the medium. 
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I. INTRODUCTION 

The correct definition of electromagnetic flux has long 
been controversial, with the main competition being be- 
tween the Abraham ^ x [H, Q and Minkowski z3 x i? Q 
forms. Pfeifer et al. [3| gave an excellent discussion and 
historical review of the situation, with an analysis based 
on energy-momentum tensors. In contrast, the Poynting 
theorem Q and Poynting vector S = E x H, in concert 
with the electromagnetic energy density, lead us to an 
energy continuity equation. This equation is easily inter- 
preted when considering only fields in the vacuum or in 
nondispersive linear media [6] . Further alternatives exist 
where the Poynting vector is generalized to include ex- 
tra terms and so generate other equally valid flux vectors 
and energy densities 0, H, S Ell ! however all these 
were based on E x H. 

However, outside the context of energy-momentum 
tensor definitions, but nevertheless common in electro- 
magnetic usage, is the E x B form [l^, [l^. Curiously, 
comparison or contrast of the Abraham/ Minkowski 
forms with E x B is. hard to find - e.g. [1] does not 
remark on the different origins of E x B at all. Further, 
even though (e.g.) both E x H and E x B appear in 
this is outside the context of magnetic media, so that 
H = B / fiQ in any case. 

In this paper we address the question: if we construct 
electromagnetic flux vectors ("Poynting vectors") using 
the cross product of two fields, what do the results look 
like, and what might we apply them to? How should 
non-trivial polarization and magnetization effects be in- 
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terpreted? Also, does the alternative D x H combina- 
tion of field vectors give interesting results? We allow for 
reasonably general propagation media, with potentially 
dispersive and nonlinear properties affecting both elec- 
tric and magnetic fields. Such complications mean that 
we do not derive continuity equations with a perfectly 
balanced fiux and energy density, but they also include 
extra "residual" terms: e.g. the standard E x B deriva- 
tion produces a residual "work done" term E ■ J, where J 
is the total current [l^[l^ . We show that in each case the 
residuals contain either temporal or spatial derivatives, 
and that they can be interpreted in terms of currents; we 
further show how they affect the propagation of the flux 
vector. 

Since typical response models for a propagation 
medium, especially in the nonlinear case, are likely to 
be non-covariant, we tie our description to the medium 
rest frame without further loss of generality. Although 
this is a restriction we might prefer to avoid, as done 
for negative refraction by McCall [l^ , it usually has few 
consequences. Further, since our interest is primarily on 
propagating fields in non-trivial media, we leave consid- 
eration of interfaces, as well as surface and volume inte- 
grals, for later work; likewise we do not address the uses 
of the vector potential A (or its dual [13, Ell) this 
context. 

The paper is organized as follows: in section |TT] we in- 
troduce Maxwell's equations in the forms in which they 
are used in this paper, in section [TTTl we briefly remark on 
continuity equations, and in section llVl we define our four 
electromagnetic fiux vectors and their associated conti- 
nuity equations. In section |V] we show how the residual 
terms modify the underlying propagation of the fiux vec- 
tors, and we discuss all the results in section IVTl Finally, 
in section fVIII we present our conclusions. 

Our presentation is pedagogic in that most of the dis- 
cussion is at the level of an undergraduate course in elec- 
tromagnetism. The material presented could form ex- 
tension work at the point when Poynting's theorem is 
introduced. Those aspects concerning the propagation 
of the fiux vector will be of most interest to specialists. 
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II. MAXWELL'S EQUATIONS 

Maxwell's equations for the electric field E and mag- 
netic field i? in a medium are 



1 1 1 

\7 -E = —pb H pf = —p 

Co Co eo 

V-B = 

V xE = -dtB 

V xB poJb + poJf + poeodtE, 



(1) 

(2) 
(3) 
(4) 



where (pb, Jb) and iPf,Jf) are respectively the bound and 
free (charge, current) densities. As an alternative, we can 
define an electric polarization P and magnetization M, 
and 



Jb = Jp + Jm = dtP + \/ xM 
Pb = -^-P 
D = eoE + P 

H = —B-M. 
Mo 

These allow us to rewrite Maxwell's equations as 

V-S = 

\IxE= -podt (i? + Af) 

V xH = Jf + dtD = Jf + dt (eoE + Py 



(5) 
(6) 
(7) 

(8) 



(9) 
(10) 

(11) 
(12) 



We can even rewrite eqn. pip in the unconventional 
form 



V xZ? = -eoMo^t i^H + AIj + \7 xP 
= -eoHodtH - coPoKb, 
where we have defined 

1 



Kb^Kp + Km 
ab^-V -M. 



■ V xP + atM, 



(13) 
(14) 

(15) 
(16) 



This Kb appears in the same place as a monopole cur- 
rent would if such were allowed; Ub is the bound mag- 
netic pole density. Note that Kb and Ub are merely a 
way of representing the (local) material response; we are 
not claiming that some process actually generates true 
magnetic monopoles inside the material Strictly 
speaking, this is also true of the bound electric charge 
and its currents - they are a mechanism used solely to 



represent the behaviour of the medium. Further, and just 
as for the ficticious bound electric charge density, the fic- 
ticious bound monopole density necessarily integrates to 
zero over all space. Thus the material response could, in 
principle, be re-represented as magnetic dipoles instead 
of monopoles. Note that in using this effective monopole 
current, we are not going as far as Carpenter who 
posits a complementary universe dominated by magnetic 
monopoles and with no charge in order to clarify some 
conceptual difficulties. If we were to include free mag- 
netic monopoles and free magnetic monopole currents, 
then the continuity equations given below would exhibit 
a great deal more symmetry on exchange of electric and 
magnetic effects. 

At this point is is worth noting that the equations 
above represent the effect of electric polarization in one 
of two ways: either as a current of bound charges (i.e. 
Jp ~ dtP), or as a result of bound monopole current 
loops (i.e. Kp = V xP). They also represent the effect 
of magnetization similarly: either as a current of bound 
monoples (i.e. Km = dtM), or as a result of bound 
electric current loops (i.e. Jm = V xM). Our aversion 
to free magnetic monopoles, and the widespread accep- 
tance of free electric charges may bias many readers to- 
ward an electric current picture involving Jp and Jm - 
but since these comprise bound charges which are merely 
a convenient fiction, there is no physical reason not to 
consider using bound monopoles, or even a mix of the 
two, if we have sufficient reason. Indeed, if we (micro- 
scopically) model the magnetization as arising from some 
field-induced or environmental distortion of a unit cell or 
molecule, there is little a priori reason not to model the 
magnetization as induced magnetic dipoles ~ we need not 
take the extra step of assuming the dipoles arise from 
some induced current loop. 



III. CONTINUITY EQUATIONS 

One of the major uses of flux (Poynting) vectors is 
in energy continuity equations, where we can examine 
the balance between flux and local storage of the energy 
in a medium. In electromagnetism, the flux vector that 
is usually chosen is the Abraham form of the Poynting 
vector E X H, although in some contexts the Minkowski 
form D X B is chosen. Some authors prefer an E x B 
form for the Poynting vector (e.g. the recent [l3|); but 
only in media with a magnetic response does this differ 
from the Abraham form in anything but scaling. 

Our starting point is just a cross-product of two se- 
lected fields, one electric (E or D) and one magnetic {H 
or B). In concert with Maxwell's equations, such cross- 
products result in continuity equations of the form 



V •S' = dtU + S. 



(17) 



^ Chapter 9, section 3 



Here S is an energy fiux based on our chosen pair of 
field vectors. The energy density U is a function of the 



2 



EMC0N4 



Four Poynting theorems 



|D r.Paul.Kinsler@physics.org| 
http: / / www.kinsler.org/physics/f 



fields used to construct the energy flux vector, and will in 
general contain all the terms that can be expressed as a 
simple time derivative of some function. The remaining 
term S is some residual contribution. One way of avoid- 
ing these residual terms is to use an energy flux vector of 
the Umov form, i.e. as an energy density multiplied by a 
velocity vector, and not the traditional cross-product of 
fields (see e.g. In the Umov approach, any terms 

that are not part of the defined energy density are at- 
tributed to the behaviour of the energy flux. 

We might now assume that the fields and medium are 
in some nearly steady state, where the time dependence 
of U is either zero, or its rapid oscillations average to zero; 
but that the residual terms remain significant. Here, any 
change in energy fiow S is balanced by the residual terms 
S which would most simply be a rate of energy removal 
or supply. However, it is perfectly possible for there to 
be more complicated responses, with a dependence on 
either time or space; e.g. an in-effect temporary storage 
of energy giving rise to an oscilliatory behaviour. This is 
treated more rigorously in section |Vl 



IV. FLUX VECTORS 

Here we will use the different forms of Maxwell's equa- 
tions given above to generate four different electromag- 
netic energy continuity equations. The way we generate 
these continuity equations is straightforward: we take 
our chosen flux vector as deflned by a cross product of 
an electric field and a magnetic field and take the diver- 
gence, using the standard vector identity 



X xY 



X ■ Vxr +Y ■ [WxX 



(18) 



Then, by substituting in appropriate Maxwell's equations 
into the RHS to substitute for the curl terms, we gen- 
erate continuity equations of the general form given in 
eqn. (flT)) . With this style of derivation, there are some 
subtleties regarding the role of total and self-field contri- 
butions, as has been discussed by Campoz and Jimenez 

Lastly, if we so wish, we might also apply the modi- 
fications to our chosen Sxy — X x Y flux vector that 
have been applied to the (bare) Abraham E x H form 
0, H, H, [TqI, El - e.g. adding terms of zero divergence to 
S 



equation, i.e 
V 



(^E X = -E ■ {v xH^ + H ■ {v xE^ (19) 
= ~E ■ dtD - E ■ Jf - H ■ dtB. (20) 



This is in itself a widely used expression, but we proceed 
further to get 

■ (e X ^ -E ■ (.Jf + (odtE + dtP) - H-fio (dtH + dtM^ 

(21) 

- -\dt (eoE ■ E + fioH ■ H 

- E ■ Jf - E ■ dtP - fioH ■ dtM (22) 

V -Seh = -dtUEH -E-Jf-E-dtP-^ioH- dtM, 

(23) 

where Seh = E x H and 2Ueh = ^qE ■ E + ^qH ■ H. 

Because we chose the fields E and H to generate our 
flux vector, we necessarily find that both of these lo- 
cal residual excitations depend only on the temporal re- 
sponse of the medium (i.e. they are dispersive). This 
means we can replace dtP with the ordinary electric cur- 
rent Jp, and dtM with a monopole current contribution 
Km- 



V -Seh = -dtUEH - E ■ Jf ~ E ■ Jp ^ ^loH ■ Km- 



XY- 



(24) 



We cannot somehow introduce Jm —VxM here without 
making assumptions: Jm because that relies on a pre- 
existing magnetization that varies in space. 

In optics, researchers work from a starting point that 
is primarily concerned with dispersion; even if that dis- 
persion has usually been purely dielectric in origin. It is 
natural, therefore, for the optics community to prefer the 
E X H form, because it treats material polarization and 
magnetization in a purely temporal manner - although 
we see here that we are then forced to represent any mag- 
netic response by bound monopoles. In the introduction 
we suggested that the E and H fields could be regarded as 
the bare fields - and it is this optics inspired temporally- 
centric view in which this is true. Another nice feature of 
this choice is that the LHS is a purely spatial derivative 
(being a divergence V •), whereas on the RHS both the 
non- Jj terms (i.e. energy density term and residuals) are 
temporal derivatives. 



A. Abraham E x H 

This Abraham form of electromagnetic flux (Poynting) 
vector is the most widely used of all, consisting of the 
fields E and H . Inserting E x H into the identity eq. 
(fT5)l . and using eqns. ^ and gives us a continuity 



Instantaneous responses 

For the case of a medium with an instantaneous scalar 
response, we have P = Xp^, so that E-dtP = {dtE-P)/2, 
and so the polarization P can be incorporated into the 
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energy density; the same can be done for the magnetiza- 
tion M. Hence we have 



1 



{ExH)=--d, 



eoil + Xp)E-E + ^0 (1 + Xm) H 
E ■ Jf (25) 



E- D + H ■ B 



E-Jf. 



(26) 



A consequence of this is that we can move any instan- 
taneous linear part of the medium response into the en- 
ergy density, leaving only the dispersive contributions to 
remain as residual terms. In practical terms, we can re- 
place eq s-iid /ip in eqn. (|23p and associated definitions 
by instantaneous-response parameters — eo(l + Xp) 
and iii — /io(l -l- Xm), then only the non-instantaneous 
or nonlinear responses will remain in P and M. 



B. Electric current: E x B 

This definition uses the electric field E and the mag- 
netic field B (see e.g. [H, 13), and has been recently 
favoured by some authors for use in complex media 
[isl [2^ . [23l . [2^ . However, one instance [23| has suf- 
fered significant (although initially disputed) comment 
[25I [26l.l27| as to its validity in certain situations. Leav- 
ing aside this dispute, we can still use the vector identity 
above to get a continuity equation; this will still be physi- 
cally valid, even if using it in complex media may mislead 
the unwary. 

To reiterate our earlier point - rather than attempting 
to justify a particular choice of fiux vector with respect to 
some external criteria, we simply set up the definition(s), 
and determine what can be done on that basis. The 
continuity equation for E x B is based on eqns. and 
dH), and is 



(^E X B^ 



-E ■ [V xBj + B ■ [V xEj (27) 
-E ■ (pLoJ + fioeodtE) ~ B ■ dtB (28) 



Mo 



dt[eoE-E 



1 

Mo 



\/-SEB = -dtUEB-E-J, 



B- B] - fioE ■ J 
(29) 
(30) 



where of course J = Jf + dtP + V xA/; this is just the 
total current density; also (J-qSeb — E x B and 2'\1eb — 
tQE-E^B-BliiQ. 

The simple representation of the residual terms above 
(i.e. as E ■ J) is presumably the underlying reason why 
E X B is preferred by many authors: the material re- 
sponse is represented in terms of an electric current, and 
E ■ J can be interpreted simply as the work done on 
charges. However, from a more general point of view, 
E X B is not the "correct" form of the Poynting vector. 



it just is one of many, each of which may (or may not) be 
more convenient in a particular situation. Note that here 
only the polarization residual excitation depends on the 
77 temporal response (i.e. dispersion); the magnetic part is 
purely spatial: both forms combine to generate a bound 
electric current. This is obviously an advantage to those 
who prefer to think only in terms of electric charges, 
whether real or ficticious; certainly E x B is the obvious 
choice for any microscopic model involving real charges 
(and no real monoples). However, the time response (dis- 
persion) of the magnetization is no longer explicit, as it 
was for E X H, but has become hidden inside the bound 
magnetization current Jm- 

In the previous subsection, we suggested that the E 
and H fields could be regarded as the bare fields. This 
is a point of view often taken in optics, and is consis- 
tent with the Ex H choice of Poynting vector. However, 
choosing Ex B gives us an alternative electric charge (or 
current) centric view, in which case it is the E and B 
fields which are regarded as the bare fields, indeed they 
are already widely regarded as the fundamental electro- 
magnetic fields Ha]^ and 



Instantaneous and pointlike responses 

In the case of instantaneous polarization response and 
pointlike magnetic response, we find that eqn. (j27|) can 
be directly reduced to the same result as in eqn. (j26p . 
but scaled by /io, since H has been replaced by B in 
the fiux vector. A consequence of this is that we can 
move such components of the medium response into the 
energy density, leaving only the dispersive polarization 
and spatial magnetic parts as residual terms. 



C. Magnetic current: D x H 

This alternate, and little (or never) used definition 
comprises the magnetic induction field H and the dis- 
placement field D as its basic components. Note that the 
historical review of Buchwald [28!| makes some discourag- 
ing remarks as regards choosing D and H as fundamental 
fields at the end of chapter 2 (p. 18), but makes no ref- 
erence to a flux generated from them. Alternatively, the 
projection approach taken to pulse propagation by Kole- 
sik et al. [1^, [s^l relied on the D and H flelds, although 
they did not consider media with a magnetic response. 
The D X H continuity equation is based on eqns. (|12p . 



^ Chapter 27, section 3 
^ Chapter 1, section 1 
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p4)) , and is 
V • (z3 X i?) 



-D- [\/xHj + H ■ {\/xD 
-D ■ (jf + dtD 



(31) 



DH 



- H ■ [fo^iodtH + eofiodtM - V xP 

(32) 

-jdt (J-3-D + noH-H 

-D-Jf- foUoH ■ dtM + H-VxP 

(33) 

-dtUoH - -D ■ Jf - fioH ■ dtM + -H 
eo Co 
(34) 



where eoSoH = D x H and 2Udh = D ■ D/e^ + fioH ■ 
H. Note that only the magnetization residual excitation 
depends on the temporal response (i.e. dispersion); the 
polarization part is purely spatial; as a result the material 
response can be encoded solely by means of the bound 
monopole current Ki, defined in eqn. (|15p : 



V -Sdh = -dtUoH ~ -D ■ Jf - ^i^H ■ Kb, (35) 

although the free current part still mimics the E ■ J form. 

This D X H form, therefore, is the natural complement 
to the electric current based Ex B form, because the ma- 
terial response in both is completely encoded by means 
of a current: but one is a magnetic monople current K^, 
and one an ordinary electric current Jb- It is probably 
unsurprising, therefore, that D x H has been neglected, 
because it treats material polarization and magnetization 
purely as monopole currents - even though these are ficti- 
tious monopoles, bound in dipole pairs, introduced solely 
to represent (model) the material response. Indeed, the 
description here is a purely continuum one, and the use 
of a monopole current is not contingent on the existence 
of fundamental particles carrying magnetic charge. 

The E X H Poyting vector gave us a time-centric view- 
point, and was compatible with regarding E and H as the 
bare fields, similarly the electric-current centred view- 
point oi E X B led us to insist that E and B are the 
bare fields. Here we have seen that a monopole-current 
centred viewpoint might encourage the idea that D and 
H are the bare electromagnetic fields! 



Pointlike and instantaneous responses 



the energy density, leaving only the dispersive magnetic 
and spatial polarization parts as residual terms. 



D. Minkowski D x B 

This uses the usual definition of the Minkowski Poynt- 
ing vector, which contains the displacement and magnetic 
fields D and B. Its continuity equation is based on eqns. 
(g]) and m]), along with eqns. ^ and (O, and is 



y -{Dx B 
V xP, 



-D- [V xB^j + B-[y xD 
-D ■ fiol^Jf + S/ xM + dtD 
+ B ■ (-eodtB + \7 xP 



1 



1 



V -Sdb = 



(36) 



(37) 



2 V eo A^o 
- noD ■ l^Jf + S/ xM^ + B - "^ xP 

(38) 

-dtUoB -—D-(jf + S/ xm) + —B • V xP, 



dt —D ■ D + —B ■ B 



(39) 



where coHoSdb — DxB and 2Udb = D-D/eo+B-B/^o- 
Note that neither residual excitation depends on the 
temporal response (i.e. dispersion); the polarization and 
magnetization parts are purely spatial. This means that 
we might replace V xM with the electric current Jm, 
and V xP with the monopole current Kp/eofiQ: 



V -Sdb = -dtUoB --D-[Jf + Jm 



B-Kp. (40) 



This form, therefore, treats material polarization and 
magnetization in a purely spatial manner; as such it pro- 
motes a picture wherein it is the D and B fields which 
look like the bare fields. 



Pointlike responses 

In the case of pointlike magnetic and polarization re- 
sponses, we find that eqn. (j36p can be directly reduced 
to the same result as in eqn. (1261) . but scaled by eo/io, 
since E and H have been replaced by D and B. A con- 
sequence of this is that we can move such components 
of the medium response into the energy density, leaving 
simplified residual terms. 



In the case of pointlike polarization response and in- 
stantaneous magnetic response, we find that eqn. f3T|) 
can be directly reduced to the same result as in eqn. 
(PS)) . but scaled by eq. A consequence of this is that we 
can move such components of the medium response into 



V. PROPAGATION OF FLUX 

We can now show how these residual terms affect a 
propagating wave, by deriving propagation equations for 
the Poynting vectors themselves. To do this we use the 
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concept of directional electromagnetic fields [3l|, |33, l33| , 
and as a result do not need to resort to restrictive approx- 
imations, such as e.g. assuming plane wave or harmonic 
fields. Notably, we will not need to resort to ad-hoc time 
averaging of fast-oscillating terms, as done in this kind 
of context by e.g. Markel and others [ISjIsSl; neither do 
we need to introduce pulse envelopes, co-moving frames, 
or make smoothness assumptions ^35, 36, 37, 38]. 

First we note that each of our electromagnetic conti- 
nuity eqns. §^,ha.s the general form 



-cV-{X xY] =^dt 



X-X + Y -Y 



R. 



(41) 



Now assume transverse fields propagating in the direction 
given by a unit vector u, so that X ■ u = Y ■ u = Q. 
This means we can construct directional fields [3l|, [s^l 
by defining 



= X^uxY- 



(42) 



where e.g. for plane polarized fields we might use 
Xx = y/f^Ex = Dx/^ and Yy = ^/Ji^Ry = By/^/JI^; 
however note that we need not be restricted to scal- 
ing by the vacuum values of eg and /iq. Thus since 
X xY = m(|G+|2 - |G^P)/4, we can write 







2 




2 








2 




2" 




G+ 












G+ 


+ 


G- 





-c V -u 



which is easily rearranged to 



-i?, 
(43) 



2 2 

[dt -cV-u] G+ + [dt + cV-u\ G" = -AR. (44) 



This contains two counter-propagating components, each 
evolved by its own wave operator 9( =F c V -u, with the 
residual terms remaining on the RHS. If we are only in- 
terested in unidirectional propagation we can set G^ — Q, 
and with u\\ z, we get 



[dt - cd,] G+ = -4R 



(45) 



This is a simple first order wave equation for the intensity, 
and we can easily see that loss-like residual terms R will 
cause that intensity to diminish; similarly those depen- 
dent on the past can cause dispersion. However, in the 
same way as with standard directional fields [3l| or fac- 
torization approaches [H, , we require that the resid- 
ual terms only have a small effect over the scale of one 
wavelength in order for G^ to stay negligible [U \^ [4^ . 

Note that a considerable amount of detail is hidden in 
the total residual term R; the various residual compo- 
nents are summarized in table U for each combination of 
field vectors. If we wanted to solve eqn. ((45|) (or perhaps 
even (|44|)). we would need to rewrite those residual com- 
ponents in terms of G^ ; the approximation G~ = helps 
in this respect since it allows the field choice X to be ex- 
pressed in terms of Y. In any case, we can see that the 
components of R will act as source terms that drive and 
modify the otherwise simple linear wave propagation. 



TABLE I: Summary of flux vectors and their corresponding energy densities, along with the residual terms and their corre- 
sponding bound currents. 



Flux 


Energy density 


Temporal 


Spatial 


Currents 


Ex H 


eoE-E + noH ■ H 


E ■ dtp + H ■ dtM 





Jp,Km 


Ex B 


eoE ■ E + fio^B ■ B 


E-dtP 


+E ■ V xM 


Jb 


DxH 




H-dtM 


-H-VxP 


Kb 


D X B 


e^^D ■ D + fig^B ■ B 





~B - V xP + D - V xM 


Kp, Jm 



VI. DISCUSSION 

Here we compare and contrast the residual terms, 
which neither appear in the form of an energy flux nor 
an energy density; they are summarized in table IH The 
Abraham E x H form has residuals that are purely 
dispersive, i.e. depend on the temporal derivatives of 
the medium responses via P and M. In contrast, the 
Minkowski D x B form has residuals that are purely spa- 



tial i.e. depend on the curl of the medium responses. The 
lesser used forms E x B and DxH have mixed resid- 
uals - the E X B residuals conveniently match the form 
of the usual definition of the total current; whereas the 
DxH residuals match the less conventional picture of a 
material response described by bound monopoles. 

Let us now consider the effects of these residual terms 
in all four cases, where we assume that Jf = 0. 

For E X H, the continuity eqn. shows that the 
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flux vs energy density balance is modified by two residual 
terms, i.e. the local dielectric excitation E-dtP, and local 
magnetic excitation H -dtM . Both of these are temporal, 
so if we choose to propagate forward in space, which is 
a common choice when considering propagation in optics 
(see e.g. [H, HI]), during each step forward in space, the 
fields E(t), H{t) will be modified according to their time 
derivatives. 

For E X B, the continuity eqn. (|^^ shows that the 
flux vs energy density balance is modified by two residual 
terms, one temporal E ■ dtP, and one spatial E ■ V xM; 
in concert they represent the local electrical work done 
on the charges comprising the current J. Unlike in the 
optics E X H picture, here a propagation step forward in 
space would be complicated by the spatial term; likewise 
the alternate choice of a time propagation step would be 
complicated by the temporal term. 

For D X H, the continuity eqn. (f33|) shows that the 
flux vs energy density balance is modified by two residual 
terms, one spatial H -W xP, and one temporal —H -dtM; 
in concert they represent the local magnetic work done 
on the monople current K . Just as for the Ex B picture, 
this leads to complicated methodology for propagation of 
the fields. 

For D X B, the continuity eqn. ([55]) shows that the 
flux vs energy density balance is modified by two residual 
terms, both spatial: —B • V xP, and D ■ V xAI. Just as 
for the ExH picture, both residuals are of the same type 
(albiet spatial, not temporal), so that choosing to prop- 
agate in time gives us the simple case where the fields 
D{r),H{f) will be modified according to their spatial 
derivatives. Despite this, propagation treatments that 
evolve the fields forward in time (e.g. FDTD 44], or 
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[13, SEl) generally tend to persist with the use of E and 
H. 



VII. CONCLUSION 

We have summarized four distinct electromagnetic 
continuity equations, each being derived from (and be- 
ing consistent with) Maxwell's equations and the stan- 
dard constitutive relations. Each handles the electric or 
magnetic response of the medium in a slightly different 
way, so the most appropriate form needs to be chosen 
according to the system under study. 

Not only have we presented the different interpre- 
tations motivated by each of the Abraham {E x H), 
Minkowski (DxB), and standard electromagnetic ( "elec- 
tric current" ) ExB forms, but we also consider the alter- 
native ("magnetic current") Dx H form as well. Finally, 
we showed how a directional fields approach [3l| could be 
used to generate a propagation equation for each flux vec- 
tor. This propagation equation makes a relatively small 
number of assumptions, and so not only enables further 
interpretation of each flux and its interaction with the 
medium, but has potential applications in its own right 
- even in materials with a complex response (46l | , and for 
wideband or ultrafast optical pulses [47| . 



Acknowledgments 

The authors acknowledge financial support from the 
Engineering and Physical Sciences Research Council 
(EP/E031463/1). 



[1] M. Abraham, Rend. Circ. Mat. Palermo 28, 1 (1909). 
[2] M. Abraham, Rend. Circ. Mat. Palermo 30, 33 (1909). 
[3] H. Minkowski, Math. Ann. 68, 472 (1910), reprint of 

Minkowski, H., 1908, Nachr. Ges. Wiss. Goettingen, 

Math.- Phys. Kl. 1908, 53. 
[4] R. N. C. Pfeifer, T. A. Nieminen, N. R. 

Heckenberg, and H. Rubinsztein-Dunlop, Rev. 

Mod. Phys. 79, 1197 (2007), 710.0461, URL 

|http : //li nk . aps . org/doi/10 . 1103/RevHodPhys .79. 1197! 
[5] J. H. Poynting, Proceedings of the Royal 

Society of London 38, 168 (1884), URL 

http : //rspl . royalsocietypublishing. org/ content /38/235- 236(^168 . full .pdf +html 
[6] L. D. Landau and E. M. Lifshitz, Electrodynamics of 

Continuous Media (Pergamon, Oxford and New York, 

1984). 

[7] J. Slepian, J. Appl. Phys. 13, 512 (1942), URL 

http://link.aip.Org/link/7JAPIAU/13/512/l 
[8] C.'S.' Lai, Am. J. P hys. 49, 841 (1981), URL 

|http : //dx ■ doi . org/ lOTl l 19/1 . 12719 
[9] P. C. Peters, Am. J. Phys. 50, 1165 (1982) 



|http://dx. doi ■org/10.1119/1 .129031 

[11] D. H. Kobe, Am. J. Phys. 50, 11 62 (1982), URL 
| http://dx. doi. org/10. 1119/1. 12901 I 

[12] R. P. Feynman, R. B. Leighton, and M. Sands, Feynman 
Lectures on Physics, The, vol. 2 (Addison- Wesley, 1964), 
14th printing. 

[13] D. B. Melrose and R. C. McPhedran, Electromagnetic 
Processes in Dispersive Media (Cambridge University 
Press, Cambridge, 1991), ISBN 0-521-41025-8. 
[14] J. R. Reitz, F. J. Milford, and R. W. Christy, Foundations 
of Electromagnetic Theory (Addison- Wesley, 1980), 3rd 



[15] F. Richter, M. Florian, and K. Henneberger, Eu- 

rophys. Lett. 13, 117 (2008), 07 10.0515, URL 

| http: // www. iop.org/EJ/abst ract/-li nk=8752695/029 5-5075/81/e 
[16] M. W! McCaU, Metamaterials 2, 92 (2008), URL 

http://dx.doi.org/10. 1016/ j .metmat . 2008 . 05 . 001 
[17] M. Hillery and L. D. Mlodinow, 
Phys. Rev. A 30, 1860 (1984), URL 



. ,, URL 

[http : //dx . doi . org/10 .1119/1 . 13024 
[10] R. H. Romer, Am. J. Phys. 50, 1166 (1982), URL 



ht tp : //link . aps . org/abstract/PR A/v30/pl860 
[18] S. Datta, Eur. J. Phys. 5, 243 (1984), 



URL 



|http://www.iop.org/EJ/abstract/0 143-0 807/5/4/010/| 



7 



EMC0N4 



Four Poynting theorems 



|D r.Paul.Kinsler@physics.org| 
http: / /www.kinsler.org/physics/f 



[19] 

[20] 

[21] 

[22] 

[23] 

[24] 

[25] 
[26] 

[27] 

[28; 
[29; 

[3o; 

[3i; 



C. J. Carpenter, lEE Proceedings - Science, Mea- http://link.aps.org/abstract/PRA/v72/e063802! 
surement and Technology 146, 73 (1999), URL [3 4] A. Favaro, P. Kinsler, and M. W. McCall, Opt. Express, 
[http : //ieeexplore . ieee . org/search/wrapper . jsp?aj:n-umber =7665S^pear (2009):comment Comment on Markel, Opt. 



M. E. Crenshaw and N. Akozbek 
Phys. Rev. E 73, 056613 (2006), URL 
|http : //link . aps . org/doi/lO . 1 103/PhysRevE . 73 . 056613 
I. Campos and J. L. Jimenez, Eur. 
J. Phys. 13, 117 (1992), U RL 

[http : //www . iop . or g/EJ/ref s/014 3-0807/ 13/3/003l 
N. Obukhov and F. W. Hehl, Phys. Lett. A 



Y 



URL 



Express 16, 19152 (2008). 
[35] R. W. Boyd, Nonlinear Optics (Academic Press, New 

York, 2003), 2nd ed., 1st ed. 1994. 
[36] T. Brabec and F. Krausz, Phys. 
Rev. Lett. 78, 3282 (1997), URL 
|http : // link . ap s . org/a'bstract/PRL/v78/p3282 
[37] M. Geissler, G. Tempea, A. Scrinzi, 
M. Schniirer, F. Krausz, and T. Brabec, 
Phys. Rev. Lett. 83, 2930 (1999), URL 
http : //link . aps . org/abstract/PRL/v83/p2930 
[38] P. Kinsler and G. H. C. New, Phys. Rev. A 

67, 023813 (2003), ph ysics/0212016vl, URL 

|http : // link . aps . org/abstract/PRA/v67/e023813 ' 
[39] P. Kinsler (2007), "Pulse propagation meth- 
http : //www. opticsinfobase . org/ abstract . cfm?URI=oe-16-23-d(M52 in nonlinear optics", 0707.0982, URL 
L. P. Pitaevsku, Phys. Rev. A 73, 047801 (2006), URL http://arxiv.org/abs/0707.0982 

|http : //link. aps . org/doi/10 ■ 110 3/PhysRevA . 73 . 047801 [40] G. Genty, P. Kinsler, B. Kibler, and J. M. 
C. Raabe and D.-G. Welsch, Phys. Dudley, Opt. Express 15, 5382 (2007), URL 

Rev. A 73, 047802 (2006), URL |http : //www, opticsexpress . org/abstract . of m?id=132608| 

|http : //link . aps . org/d oi/10 . 11 03/PhysRevA . 73 . 047802 ' [41] P. Kinsler, J. " 



311, 277 (2003), (from Brevik-E-2008pra) 
[dx.doi. org /10 . 1016/S0375-9601 (03) 00503-6 

C. Raabe and D.-G. Welsch, Phys. 

Rev. A 71, 013814 (2005), URL 
|http : //link . aps . org/doi/10 . 1103/PhysRevA . 71 . 013814] 

V. A. Markel, Opt. Express 16, 

19152 (2008), 0712.0605, URL 



Opt. Soc. Am. B 24, 



2363 



I. Brevik 
Rev. A 



and 
79, 



S. A. 

027801 



Ellirigsen, 
(2009), 



Phys. 
URL 



[32; 
[33; 



http : //link. aps . org/doi/10 . 1103/PhysRevA . 79 . 027801 

J. Z. Buchwald, From Maxwell to Microphysics (Univer- 
sity of Chicago, 1985), ISBN 0-226-07883-3. 
M. Kolesik, J. V. Moloney, and M. Mlejnek, 
Phys. Rev. Lett. 89, 283902 (2002), URL 
http : //link. aps . org/abstract/PRL/v89/p283902 
M. Kolesik and J. V. Moloney, Phys. 
Rev. E 70, 036604 ( 2004), URL 
|http : //link . aps . org/abstract /PRE/ v70/e036604 
P. Kinsler, S. B. P. Radnor, and G. H. C. New, 
Phys. Rev. A 72, 063807 (2005), note that in this 
reference, the convolution symbol between the acPr 
coefEcents and the terms in square brackets 

in the frequency-domain propagation equations was 
inadvertently omitted, see physics/0611215vl, URL 
|http : //link . aps . org/abstract /PRA/v72/e063807 
P. Kinsler (2006), "Theory of directional pulse prop- 
agation: detailed calculations", physics/0611216, URL 
http : //arxiv. org/abs/physics/0611216 
Y. Mizuta, M. Nagasawa, M. Ohtani, and M. Ya- 
mashita, Phys. Rev. A 72, 063802 (2005), URL 



(2007), the arXiv:0707.0986 version contains 
an additional appendix., 0707.0986v2, URL 
http: //josab. osa. org/ abstract . cf m?id=140996 
[42] P. Kinsler (2008), "Transverse hmits on the uni- 
directional pulse propagation approximation", 
0810.5701, URL http://arxiv.org/abs/0810.5701 
[43] P. Kinsler (2008), "Optical pulse propagation 
with minimal approximations", 0810.5689, URL 
|http : //arxiv . org/abs/0810 . 5689 ; 
[44] K. S. Yee, IEEE ' Trans. Anten- 

nas Propagat. 14, 302 (1966), URL 

|http: / /ieeexplore . ieee . org/xpl/f reeabs_all . jsp?arnvunber=113i 
[45] M. Scalora and M. K Crenshaw^ 
Opt. Commun. 108, 191 (1994), URL 
[http:// dx . doi . org/10 . 1016/0030-4018 (94) 90647-5 
[46] Focus Issue: Negative refraction and metamate- 

rials. Opt. Express 11, 639-760 (2003), URL 

|http: //www\^pticsinf obase . org/oe/issue . cfm?voluiiie=ll&issue= 
[47] T. Brabec and R Krausz^ Rev. 
Mod. Phys. 72, 545 (2000), URL 
|http : //prolaLTaps . org/abstract /RMP/v72/i2/p545-lj 



8 



